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This study examines the impact of parameterization of two variables, light use efficiency (LUE) and the 
fraction of absorbed photosynthetically active radiation (fPAR or fAPAR), on gross primary production 
(GPP) modeling. Carbon sequestration by terrestrial plants is a key factor to a comprehensive under- 
standing of the carbon budget at global scale. In this context, accurate measurements and estimates of 
GPP will allow us to achieve improved carbon monitoring and to quantitatively assess impacts from cli- 
mate changes and human activities. Spaceborne remote sensing observations can provide a variety of 
land surface parameterizations for modeling photosynthetic activities at various spatial and temporal 
scales. This study utilizes a simple GPP model based on LUE concept and different land surface param- 
eterizations to evaluate the model and monitor GPP. Two maize-soybean rotation fields in Nebraska, 
USA and the Bartlett Experimental Forest in New Hampshire, USA were selected for study. Tower-based 
eddy-covariance carbon exchange and PAR measurements were collected from the FLUXNET Synthesis 
Dataset. For the model parameterization, we utilized different values of LUE and the fPAR derived from 
various algorithms. We adapted the approach and parameters from the MODIS MODI 7 Biome Properties 
Look-Up Table (BPLUT) to derive LUE. We also used a site-specific analytic approach with tower-based 
Net Ecosystem Exchange (NEE) and PAR to estimate maximum potential LUE (LUE max ) to derive LUE. For 
the fPAR parameter, the MODIS MODI 5A2 fPAR product was used. We also utilized fAPAR ch i , a parameter 
accounting for the fAPAR linked to the chlorophyll-containing canopy fraction. fAPARchi was obtained by 
inversion of a radiative transfer model, which used the MODIS-based reflectances in bands 1-7 produced 
by Multi-Angle Implementation of Atmospheric Correction (MA1AC) algorithm. fAPAR cU exhibited sea- 
sonal dynamics more similar with the flux tower based GPP than MODI 5A2 fPAR, especially in the spring 
and fall at the agricultural sites. When using the MODIS MOD17-based parameters to estimate LUE, 
fAPARchi generated better agreements with GPP (r 2 = 0.79-0.91) than MOD15A2 fPAR (r 2 = 0.57-0.84). 
However, underestimations of GPP were also observed, especially for the crop fields. When applying the 
site-specific LUE max value to estimate in situ LUE, the magnitude of estimated GPP was closer to in situ 
GPP; this method produced a slight overestimation for the MOD15A2 fPAR at the Bartlett forest. This 
study highlights the importance of accurate land surface parameterizations to achieve reliable carbon 
monitoring capabilities from remote sensing information. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

Carbon sequestration by terrestrial plants is a key factor to 
a comprehensive understanding of carbon budget at the global 
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scale. Carbon assimilation through photosynthesis by vegetation, 
known as gross primary production (GPP) at the ecosystem scale, 
is the largest carbon exchange between the biosphere and the 
atmosphere (Beer et al., 2010). Accurate measurements and esti- 
mates of GPP will allow us to achieve improved carbon monitoring 
and to quantitatively assess impacts from climate changes and 
human activities (Grace et ah, 2007 ; Middleton et ah, 201 1 ). Remote 
sensing observations can provide a variety of land surface parame- 
terizations for modeling photosynthetic activities at various spatial 
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and temporal scales, and much effort has been put towards this 
goal over the past decades. Many remote-sensing-based indices 
and algorithms tracking leaf biochemical properties (e.g., chloro- 
phyll, water content) and canopy biophysical properties (e.g., leaf 
area index, LAI; fraction of photosynthetically active radiation 
(PAR) absorbed by vegetation canopy (fPAR or fAPAR) have been 
developed to estimate GPP with different approaches for various 
ecosystems with variable success (Cheng et al., 2007; Grace et al., 
2007; Heinsch et al., 2006; Houborget al., 2011; Peng and Gitelson, 
2012; Running et al., 2004; Sellers et al., 1997; Tucker and Sellers, 
1986; Turner et al„ 2006; Wu et al., 2011, 2009; Xiao et al., 2010, 
2004a; Zhang et al., 2009). 

The concept of the light use efficiency (LUE) model (Monteith, 
1972, 1977) has been widely used for GPP modeling. The approach 
describes carbon sequestration as a product of the capability and 
efficiency of vegetation to capture and convert solar radiation into 
biomass by terrestrial vegetation. The model is usually defined as: 

GPP = e x PARx/PAR (1) 

where e is LUE over a defined period of time. Many previous stud- 
ies have utilized the model and coupled it with remote sensing 
data to estimate GPP at various spatial scales (Cheng et al., 2009; 
Drolet et al., 2008; Hall et al., 2012; Hilker et al., 2012; Huntzinger 
et al., 2012; Middleton et al., 2009; Nichol et al., 2000; Schaefer 
et al., 2012), including the Moderate Resolution Imaging Spectrora- 
diometer (MOD1S) Terra/Aqua satellite GPP Products (Heinsch et al., 
2006; Running et al., 2004; Turner et al., 2006; Zhao et al., 2005). 
The algorithm used to deliver the MODIS MOD17A2 GPP product 
utilizes a biome-based look-up table (LUT) to determine the poten- 
tial maximum LUE (LUE max or £ max ), and assumes down-regulation 
based on ambient air temperature and vapor pressure deficit (VPD). 

The MODIS MODI 7 product utilizes the MODIS MOD15A2 fPAR 
product in the procedure. The MODIS MOD15A2 algorithm incor- 
porates atmospherically corrected MODIS surface reflectance and 
a LUT method to achieve inversion of the three dimensional radia- 
tive transfer process in vegetation canopies. A back-up method 
based on relationships between the normalized difference vege- 
tation index (NDVI) and LAI and fPAR is used when the inversion 
method fails to deliver (Knyazikhin et al., 1999, 1998; Myneni 
et al., 2002). This approach estimates the fAPAR associated with the 
entire canopy or ecosystem, comprised of both photosynthetic and 
non-photosynthetic components. The concept of retrieving fAPAR 
for the chlorophyll-containing canopy actively involved in pho- 
tosynthesis was proposed by Zhang and colleagues (Zhang et al., 
2013, 2012, 2009, 2005). This retrieval was designed to account for 
only the solar energy absorbed by chlorophyll of the foliage and 
contributing to photosynthesis. This was accomplished by inver- 
sion of a linked leaf-canopy radiative transfer model using surface 
reflectance of five or seven MODIS land bands (Zhang et al., 2013, 
2012,2009, 2005). 

The vegetation photosynthesis model (VPM) is another simple 
GPP model based on remote sensing observations and microme- 
teorological data (Kalfas et al., 2011; Wang et al., 2010; Xiao et al„ 
2004a, b; Xiao et al., 2005). The VPM model adapted the LUE concept 
and calculates GPP as the product of PAR, fAPAR of the photosyn- 
thetically active vegetation (denoted as fAPAR pav or fAPAR ch i), and 
LUE (e). The Enhanced Vegetation Index, EVI (Huete et al., 2002), 
was used as a proxy of fAPAR ch i (Xiao et al., 2004a, b). LUE (e) is 
derived as the potential maximum LUE (e max ) down-regulated by 
scalars of temperature, water, and phenology. The e max value is usu- 
ally determined either from previous studies or from an analytical 
approach based on flux tower net ecosystem exchange (NEE) and 
PAR. This model has been used and examined with various vegeta- 
tion types including deciduous and evergreen forests, crop (maize), 
grassland, and savannah woodlands (Jin et al., 2013; Kalfas et al., 
2011 ; Wang et al., 2010; Xiao et al., 2004a, b, 2005). 


The importance of accurate estimate of both fPAR and LUE for 
improved carbon monitoring capabilities has been pointed out in 
recent studies. For instance Bonan and colleagues (Bonan et al., 
2011) suggested the accurate radiative transfer parameterization 
within the canopy (e.g., fPAR) is necessary to improve output accu- 
racy of the simulated GPP from Community Land Model version 4 
(CLM4). LUE has been shown to have significant variations across 
vegetation types (Anderson et al., 2000; Gower et al., 1999). When 
applied to models, large errors were found and stem from biases 
in LUE that accumulate in the calculation of annual carbon assimi- 
lation (Lin et al., 201 1 ). Therefore, in this paper, we integrated the 
LUE concept with different approaches for determining the param- 
eters e and fPAR to estimate carbon assimilation (GPP) at three 
FLUXNET sites, two agricultural sites and one deciduous forest site. 
Our objectives are to validate the capability to model and moni- 
tor GPP using the simple LUE concept and to compare results from 
different parameterization schemes. 

2. Methods 

2.1. Study sites 

In this study, we investigated the capability of a simple model 
based on the LUE concept to estimate GPP at three FLUXNET 
sites using different land surface parameterizations. The first two 
sites are maize-soybean rotation crop fields located near Mead, 
Nebraska, USA. The third site is the Bartlett Experimental For- 
est in New Hampshire, USA. These sites were selected based on 
their phenological characteristics and the availability of flux data, 
micrometeorological measurements, and remote sensing observa- 
tions. Furthermore, they provided us the opportunity to examine 
the model on different vegetation types in both native and managed 
ecosystems. 

2.1.1. Maize-soybean rotation crop fields, Nebraska (US-NE2 and 
US-NE3) 

Two of the three crop fields located at the University of Nebraska 
Agricultural Research and Development Center near Mead, NE, USA 
were selected to study. Both fields have maize ( Zea mays L.)- soy- 
bean ( Glycine max [L.] Merr.) rotations, but depend differently on 
water supply. The 52.4ha US-NE2 site (41.1649°N, 96.4701 W) is 
equipped with a center pivot irrigation system while the 65.4 ha 
US-NE3 site (41 .1 797°N, 96.4396°W) relies on rainfall (Verma et al., 
2005). The soils are deep silty clay loams. The fields were uniformly 
tilled by disking prior to 2001 and since then have been managed 
as a no-till system. Both fields are equipped with eddy covariance 
tower systems to measure energy and C0 2 fluxes (Verma et al., 
2005). In this study, the data from 2001 through 2004 were col- 
lected at the two sites and analyzed. 

2.1.2. Bartlett experimental forest, New Hampshire (US-BAR) 

The Bartlett Experimental Forest (44.0646°N, 71.2881°W) is 
located within the White Mountains National Forest in northern 
central New Hampshire, USA (Jenkins et al., 2007). The 1050 ha for- 
est is managed by the USDA Forest Service Northeastern Research 
Station. The forest experiences a humid continental climate with 
short cool summers and long cold winters (Jenkins et al., 2007). 
A 26.5 m high tower was installed for continuous eddy covari- 
ance fluxes and meteorological measurements in a low-elevation 
northern hardwood stand, where the average canopy height was 
approximately 20-22 m. In the tower footprint, the forest is pre- 
dominantly classified into red maple (Acer rubrum ), sugar maple 
(Acer saccharum), and American beech (Fagus grandifolia), for which 
LAI in the vicinity was reported to be 3.6 to 4.5 (Jenkins et al., 
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2007). Available data from 2004 through 2005 were collected for 
this study. 

2.2. Flux data collection 

Flux and meteorological variables were collected from the 
FLUXNET Synthesis Dataset (http://www.fluxdata.org) for the 
study sites. The Synthesis Dataset provides a collaboration por- 
tal for flux and meteorological data collected worldwide (Agarwal 
et al., 2010). The data are quality-checked and gap-filled based on 
guidelines and previous studies (Agarwal et al., 2010; Falge et al., 
2001 ; Munger and Loescher, 2006). NEE was processed following 
the FLUXNET procedure, from which GPP values were computed 
(Agarwal et al., 2010). For the maize-soybean rotation crop fields 
(US-NE2, US-NE3), daily observations, including air temperature, 
vapor pressure deficit (VPD), GPP, and PAR, were collected from 
2001 through 2004. For the Bartlett forest (US-BAR), daily mea- 
surements from 2004 and 2005 were obtained. 

2.3. MODIS data processing 

The MODIS data were processed using the Multi-angle 
implementation of atmospheric correction (MAIAC) algorithm 
(Lyapustin et al., 201 la, b, 2012). MAIAC is an advanced algo- 
rithm which uses time series analysis and a combination of 
pixel-based and image-based processing to improve the accu- 
racy of cloud/snow detection, aerosol retrievals, and atmospheric 
correction based on the BRDF model of surface. For this work, 
we used MODIS Collection 6 calibrated and geolocated (level 
IB) data. The derived surface reflectance in MODIS bands 1 
(620-670 nm), 2 (841 -876 nm), 3 (459-479 nm), 4 (545-565 nm), 
5 (1230-1250 nm), 6 (1628-1652 nm), and 7 (2105-2155 nm) were 
then utilized to calculate fAPAR ch i, as shown below, for the study 
sites. 

2.4. Model parameterization 

For the GPP model [ 1 ] utilized in this study, the PAR values were 
obtained from the FLUXNET dataset. Different approaches were 
used to determine the remaining two parameters.e and fPAR. 

2.4.1. Light use efficiency 

The MODIS MODI 7 GPP products (Running et al., 2004) come 
from an algorithm that describes the parameter e as the down- 
regulated value of the maximum radiation conversion efficiency 
(£max). The down-regulation is driven by sub-optimal environmen- 
tal conditions determined by air temperature and VPD: 

£ = £max x T min scalar X VPD scalar (2) 

where Tmin sca | ar and VPD sca [ ar are dimensionless multipliers 
describing relative decreases of e max and photosynthesis due to 
cold temperature and high VPD. These scalars are defined as simple 
linear ramp functions using prescribed minimum and maximum 
environmental constraints (Tmin mirl and rmin max , VPD min and 
VPD max ) based on broad global biome classes. Detailed descrip- 
tions can be found in previous publications (Fleinsch et al., 2003, 
2006). We adapted the Biome Properties Look-Up Table (BPLUT) 
used to derive the MODI 7 product. The values ofe max and parame- 
ters needed to calculate Tmin sca i ar and VPD sca ] ar are listed in Table 1 . 

In addition to the BPLUT values used in the MODIS MODI 7 prod- 
uct (Table 1 ), we also utilized a site-specific analysis between NEE 
and PAR to determine £ max . The analysis is usually done by per- 
forming a linear or nonlinear fitting to half-hourly or hourly data 
collected at flux tower sites (Ruimy et al., 1995; Wofsy et al., 1993). 
The approach has been adapted in many previous studies to deter- 
mine the e max value for different sites (Wang et al., 2010; Xiao, 


Table 1 

Biome properties look-up table (BPLUT) describing parameters used in GPP model. 


Site 

US-NE2 and US-NE3 

US-BAR 

Vegetation type 

Crop 

Deciduous broadleaf forest 

£max (kgC MJ-1) 

0.000680 

0.001044 

Tmin mi „ (°C) 

-8.00 

-8.00 

Tmin max (“C) 

12.02 

7.94 

VPD min (Pa) 

650 

650 

VPD m ax (Pa) 

4100 

2500 


2006 ; Xiao et al., 2004a, 2005 ). Following the description and exam- 
ples established in previous studies (Wang et al., 2010; Xiao et al., 
2004a, 2005), we utilized the hourly PAR and flux data collected 
from the FLUXNET Synthesis Dataset and applied a nonlinear hyper- 
bolic model to estimate e max for each of the study sites. The derived 
£max values were then used as input in Eq. (2) to calculate e, in addi- 
tion to the values computed based on e ma x adapted from MODIS 
MODI 7 algorithm (Table 1). 

2.4.2. Fraction of absorbed PAR 

We utilized two different values for the fPAR parameter. One 
is from the MOD15A2 LAI /fPAR product. The MOD15A2 LAI/fPAR 
product was downloaded from Oak Ridge National Laboratory 
(ORNL) Distributed Active Archive Center (DAAC) for Biogeochemi- 
cal Dynamics (http://daac.ornl.gov/MODIS/). The MODIS MODI 5A2 
fPAR product represents the fraction of PAR absorbed by the whole 
canopy. 

The concept of fAPAR associated with the foliage chlorophyll 
content (fAPAR ch |) accounts for the PAR absorbed by chlorophyll 
throughout a canopy (APAR ch i) used in the photosynthesis pro- 
cess, as compared to the PAR absorbed by the whole canopy (e.g., 
MOD15A2 fPAR). The fAPAR ch i concept partitions a vegetation 
canopy into leaf and non-leaf (referred to as stem) components. The 
leaf can be further partitioned into chlorophyll, non-photosynthetic 
pigments (referred to as brown pigment) and dry matter fractions. 
The PROSPECT (Baret and Fourty, 1997; Jacquemoud and Baret, 
1990) leaf radiative transfer model (RTM) and the SAIL2 (Braswell 
et al., 1996) canopy RTM were coupled to account for the differ- 
ent foliage fractions and full canopies. The fAPAR values for the 
non-chlorophyll foliage sector (fAPARnon.chi) and for total foliage 
(fAPARf 0 | iage ) can be computed as: 

/APAR non _ ch j = /APAR| )10wn _ pigment +/APAR drymatter (3) 

/APAR foliage =/APAR chl +/APAR non . chl (4) 

Each MODIS surface reflectance observation (p obs ) contained 
some noise, and was associated with a specific geometries for view 
zenith angle ( 8 V ), relative view azimuth angle (</>), and solar zenith 
angle ( 6 S ), all expressed in degrees. Therefore, we treated each 
reflectance observation as a sample of the following distribution: 

P ~ {Pobs(^, 0v (1 + 3N(0, 1)) , 6s (1 + 3N(0, 1 )) , </>(l + 3N(0, 1)))) 

(1 +0.05N(0, 1)) (5) 

where N(0,1) was the normal distribution with a mean of zero 
and SD = 1. We used multiple samples from this distribution in Eq. 
(5). The Markov Chain Monte Carlo (MCMC) procedure was imple- 
mented for the inversion of the models to calculate the parameter 
fAPAR ch |. A detailed description of the algorithm can be found in 
previous publications (Zhang et al„ 2005, 2009, 2012, 2013). In 
earlier works, the fAPAR dl | variable was derived using the MODIS 
MOD09 surface reflectance product and the Earth Observing-1 (EO- 
1 ) Hyperion imagery and was shown to be a useful parameter that 
improved carbon monitoring (Zhang and Middleton, 2010; Zhang 
etal., 2005, 2006, 2009, 2012, 2013). In this study, MAIAC-processed 
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MOD1S reflectance band 1-7 was used to calculate fAPAR ch i for the 
study sites. The fAPAR ch | parameter was then implemented in the 
LUE model to estimate GPP. 

In summary, we utilized two different values of e ma x values with 
the MOD1S MODI 7 algorithm, Eq. (2), to calculate s. We linked 
the e values to two different fPAR values, MODIS MOD15A2 fPAR 
and fAPAR ch i, to estimate GPP at eight-day temporal scale for all 
of the study sites. Four combinations of the different land surface 
parameterizations were examined and are summarized in Table 2. 


Table 2 

Combinations of LUE (s) and fPAR parameterization examined in this study. Note, 
values were used in the MODIS MOD17 like algorithm (Eq. (2)) to determine e 
for the GPP model. 



£max 

fPAR 

Case I 

MODIS MODI 7 BPLUT 

MOD 15A2 fPAR 

Case II 

MODIS MODI 7 BPLUT 

fAPARth, 

Case III 

Analytic approach, site specific 

MOD 15A2 fPAR 

Case IV 

Analytic approach, site specific 

fAPAR^ 



500 

450 

400 

350 ^ 
oo 

300 > 
250 ^ 
200 | 
150 £ 
100 
50 
0 



600 

500 

400 £ 
E 

300 ^ 
_E 
cc 

200 < 
a. 

100 

0 



450 

400 

350 

300 | 

250 JE 
"o 

200 _£ 
C E 

150 < 

Q- 

100 

50 

0 


Fig. 1. Seasonal dynamics of GPP (black solid line) and PAR (grey solid line) acquired from the Fluxnet Synthesis Dataset for (a) an irrigated crop site, US-NE2; (b) a rainfed 
crop site, US-NE3 in 2001 through 2004 and (c) a deciduous forest, US-BAR in 2004 and 2005. 
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Fig. 2. Seasonal dynamics of GPP acquired from the Fluxnet Synthesis Dataset (black solid line), the MODIS MOD15A2 fPAR product (■), and fAPAR ch i (A) for (a) an irrigated 
crop site, US-NE2; (b) a rainfed crop site, US-NE3 in 2001 through 2004 and (c) a deciduous forest, US-BAR in 2004 and 2005. 


3. Results 

Values of GPP (black solid line) and PAR (grey solid line) acquired 
from the FLUXNET Synthesis Dataset are shown in the format of 
eight-day values for the three study sites, rotation maize-soybean 
field (US-NE2, US-NE3) and the Bartlett Experimental Forest (US- 
BAR), in Fig. 1 a, b, c, respectively. Notable differences were observed 
for the two alternating crop types in US-NE2 and US-NE3 (Fig. la, 
b). Maize showed much higher productivity in 2001 and 2003 than 


soybean in 2002 and 2004. GPP in all three sites exhibited pro- 
nounced seasonal patterns (Fig. 1) with very low productivities 
in winter, that increased in springtime after planting (for agri- 
cultural sites) or green-up (for forest site), and with peak GPP 
observed in summer, decreasing in fall. Seasonal dynamics of the 
MODIS MOD15A2 fPAR product (■) and fAPAR dl | (A) are shown 
as eight-day average values in Fig. 2 for the three study sites. 
All three variables (GPP, MOD15A2 fPAR, and fAPAR ch |) showed 
clear differences between growing and non-growing seasons. For 
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Fig. 3. Estimated GPP using MODI 5A2 fPAR and l7\PAR d] | with e derived from the MODIS MODI 7 like algorithm for the irrigated crop site, US-NE2 (a) maize and (b) soybean. 



Fig. 4. Estimated GPP using MOD15A2 fPAR and l'APAR (h i with £ derived from the MODIS MODI 7 like algorithm for the rainfed crop site, US-NE3 (a) maize and (b) soybean. 


the two agricultural sites, US-NE2 and US-NE3, fAPAR ch | exhibited 
seasonal patterns that were more similar to measured GPP than 
those obtained with MOD15A2 fPAR. However, the MOD15A2 fPAR 
described unrealistically longer growing seasons than those shown 
by GPP or fAPAR ch i. The MOD15A2 fPAR product revealed higher 
values in spring and fall, while GPP suggested low (or no) photosyn- 
thetic activities in the field. At the Bartlett forest site (US-BAR), GPP, 
and fAPAR ch i showed a similar seasonal pattern and length of grow- 
ing season (Fig. 2c), whereas the MODI 5A2 fPAR product described 
these seasonal dynamics less accurately and exhibited much higher 
values in winter, early spring, and late fall than fAPAR ch | (Fig. 2c). 

The GPP model ( 1 ) was applied using parameterization from the 
MODIS MOD17 product algorithm (Table 1) to determine LUE (e) 
coupled with either MOD15A2 fPAR (Table 2, Case 1) or fAPAR ch i 
(Case II) to model GPP for the three study sites. Comparisons 
between in situ GPP and estimated GPP for the US-NE2 site are 
shown in Fig 3. For both maize and soybean on the irrigated US-NE2 
site, the estimates using fAPAR ch] show consistently good linear 
relationships (r 2 = 0.92- 0.94) with flux tower based GPP (Fig. 3). 
Estimated GPP using MOD15A2 product showed less satisfactory 
(r 2 = 0.61- 0.78) results, with a noticeable offset for the regres- 
sion line (Fig. 3). Nevertheless, clear GPP underestimations were 
observed with either fAPAR ch | or MOD15A2 fPAR. Results for the 
rainfed US-NE3 site are summarized in Fig. 4, where GPP estimates 
obtained using either the MOD15A2 product or the fAPAR ch i were 
comparable (r 2 =0.84- 0.87 for maize; r 2 =0.72- 0.79 for soybean). 
Stronger regressions were obtained for irrigated (Fig. 3) vs. rainfed 
(Fig. 4) fields using fAPAR chI , but the opposite for the MOD15A2 
product. Similar to the US-NE2 site, estimates using MOD15A2 
product showed a more obvious offset, and clear underestima- 
tions were observed for both fPAR products (Figs. 3 and 4). For 
the Bartlett forest, estimates using fAPAR ch i showed much stronger 


relationships to tower-based GPP (r 2 =0.91, Fig. 5) compared to 
the MOD15A2 product (r 2 = 0.57). Nevertheless, estimates using 
the MOD15A2 product exhibited slightly lower RMSE that using 
fAPAR ch |, and GPP estimates were closer to tower values than for 
fAPAR chl (Fig. 5). 

In addition to using the e ma x values adapted from the MODIS 
MODI 7 product algorithm (Table 1), we also applied an analytic 
approach based on NEE and PAR to determine the site spe- 
cific value of £ ma x- For the US-NE2 sites, 0.004209 kg CMJ -1 was 
retrieved for maize and 0.002606 kg CMJ~ ! was used for soybean. 



Fig. 5. Estimated GPP using MOD15A2 fPAR and fAPAR ch i with e derived from the 
MODIS MODI 7 like algorithm for the deciduous forest site, US-BAR. 
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in situ GPP (g C/m2/8d) in situ GPP (g C/m2/8d) 

Fig. 6. in situ GPP and estimated GPP using MOD15A2 fPAR and fAPAR d] | with e derived using site specific £ max for the irrigated crop site, US-NE2 (a) plotted against time 
from 2001 through 2004, (b) comparisons for maize, (c) comparisons for soybean. 


For the US-NE3 sites, 0.00555 kg C MJ -1 was derived for maize and 
0.00281 7 kg CMJ -1 was used for soybean. For the US-BAR forest, 
0.002006 kg CMJ -1 was used. These site specific e max values were 
used as input to the MODIS MODI 7 product algorithm, Eq. (2), to 
derive e. Flux tower based GPP and estimated GPP using the site spe- 
cific £ ma x value coupled with MODI 5A2 fPAR (Case III) and fAPAR ch i 
(Case IV) are plotted as a time series for the three study sites. Com- 
parisons are shown as regression charts in Figs. 6-8. In general, 
improvements in the magnitude of GPP estimates and reductions in 
RMSE were observed, especially for the agricultural sites. Estimated 
GPP using fAPAR chI (Case IV) was very close to flux tower based val- 
ues, for both maize and soybean (Figs. 6 and 7 a). For the Bartlett 
forest, a slight decrease in RMSE was observed for estimates using 
fAPAR ch |, but a clear increase in RMSE and GPP overestimates were 
observed for estimates using the MOD15A2 fPAR product paired 
with the site-specific £ max (Fig. 8b). 

4. Discussion 

In this study, we utilized different approaches to derive two 
important parameters needed for carbon modeling, LUE (e) and 
fPAR, with a simple GPP model. Performance of these approaches 


was assessed for simulated GPP, as compared to in situ flux tower 
GPP. 

All three study sites exhibited seasonal GPP patterns (Fig. 1 ), as 
expected. In both the maize-soybean rotation sites (US-NE2 and 
US-NE3), the C4 plant, maize, showed higher GPP values than the 
C3 plant, soybean. Furthermore, we noticed higher GPP values for 
irrigated crops (Fig. la) than rainfed crops (Fig. lb) that could be 
attributed to water availability. For the agricultural sites, little or no 
green biomass should be expected outside of the growing season. 
After planting, clear patterns of growing season were observed from 
the early vegetative stage, the mature and reproductive stage, to the 
senescence stage. Seasonal dynamics of fAPARchi closely matched 
GPP for both of the agricultural sites. The MOD15A2 fPAR product 
showed similar values to fAPAR ch | during the peak of the grow- 
ing season (approximately mid-July, Fig. 2a, b) but higher values 
for the rest of the year. The time series of MOD15A2 fPAR prod- 
uct also suggested earlier green-up and later senescence, i.e. a 
longer duration of the growing season, than determined from either 
tower observations or fAPAR ch i. At the Bartlett forest (US-BAR), 
the MOD15A2 fPAR product had higher values than fAPAR ch i not 
only outside growing season but also during the summer growing 
season (Fig. 2c). The duration of the growing season suggested by 
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Fig. 7. in situ GPP and estimated GPP using MODI 5A2 fPAR and fAPAR ch with e derived using site specific e max for the rainfed crop site, US-NE3 (a) plotted against time from 
2001 through 2004, (b) comparisons for maize, (c) comparisons for soybean. 


MOD15A2 fPAR at the US-BAR site also appeared longer compared 
with fAPAR ch i and GPP. MOD15A2 fPAR suggested a slightly earlier 
green-up and a later senescence (Fig. lc). 

For all the study sites little or no green biomass should be 
expected outside of the growing season. Flux tower based GPP also 
showed little or no photosynthetic activities between senescence 
and green-up. The fAPAR ch i exhibited seasonal variations match- 
ing the expectation. The MOD1S MOD15A2 fPAR product seemed 
to have an overestimation bias for all three sites during winter, 
early spring, and late fall. The fAPAR ch i is a parameter proposed 
here to account for the fraction of PAR absorbed by chlorophyll 
in leaves, only which will be used in the photosynthesis process. 
Theoretically, this parameter minimizes the effects caused by the 
non-green component of leaf and canopy (e.g. non-photosynthetic 
pigment and stem) as well as soil background. This might explain 
the lower values of fAPAR ch ] than MODIS fPAR product. This feature 
might also make fAPAR ch | a land surface parameter that is closer to 
the “true” physiological and phenological condition, both during 
and outside the growing season. 

The results from coupled MODIS fPAR and fAPAR ch i with the 
LUE (e) value derived from the algorithm similar to the MODIS 
MODI 7 product are shown in Figs. 3-5. Estimates using fAPAR chl 
consistently produced better agreement with in situ GPP values 


for all the three study sites. This suggests that fAPAR ch ] might 
be used as an alternative option than the standard MOD15A2 
fPAR product as input to carbon models to improve our current 
capabilities in carbon modeling and monitoring. As shown in the 
regression charts (Figs. 3-5), estimates using MOD15A2 fPAR had 
more significant offsets than fAPAR ch ]. This is most likely due to 
the higher values of MODI 5A2 fPAR exhibited outside of the grow- 
ing season. A general underestimation of modeled GPP was found 
for two agricultural sites (Figs. 3 and 4). This might suggest an 
underestimation of e values for the crops using the algorithm 
similar to the MODIS MODI 7 product. Similar underestimation 
was also found for the US-BAR sites, but was not as pronounced 
(Fig. 5). The differences in performance among sites could also 
be attributed to differences in plant functional types (C3 and C4 
crops, broadleaf trees). The model we utilized in this study took 
the ability to capture and the efficiency to use solar radiation into 
account in GPP estimates. Such plant physiological properties vary 
among different vegetation types, and might explain the differ- 
ent performance as we observed at the three sites. Similarly, a 
recent study using plant biochemical properties and environmen- 
tal conditions to estimate GPP (Wu et al., 2011) also reported that 
model calibration and performance can depend on plant functional 
types. 
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Fig. 8. in situ GPP and estimated GPP using MOD15A2 fPAR and fAPAR ct( with s derived using site specific £ max for the deciduous forest site, US-BAR (a) plotted against time 
from 2004 and 2005, (b) comparison chart. 


Using the site-specific e max values improved the accuracy of 
estimated GPP values, especially for the two agricultural sites 
(Figs. 6 and 7). As we discussed previously, the MODI 5A2 fPAR pro- 
duced GPP overestimation outside of the growing season in winter, 
early spring, and late fall for both the US-NE2 and US-NE3 sites. The 
best GPP estimates were obtained using fAPAR ch i coupled with e, 
determined by the site specific e max value and the MODIS MODI 7 
algorithm (Case IV). At the Bartlett forest, using site specific e max 
also improved GPP estimates using fAPAR ch i. However, GPP over- 
estimates resulted with the site-specific e max and with MOD15A2 
fPAR. One may notice from the regression chart (Fig. 5), that GPP 
estimates using MOD15A2 fPAR paired with MOD17 LUE (Case I) 
were close to in situ GPP values when GPP values were high, most 
likely during peak of the growing season. This suggests that the sim- 
ple model of using MOD15A2 fPAR and e values might have good 
performance for some study sites during the peak of the growing 
season. One should also note that the e max value can be obtained 
from literature review or derived using the analytical approach, 
which is affected by multiple parameters in the process including 
regression model type and the period of time of interests (Ruimy 
et al., 1995; Wang et al„ 2010; Xiao et al„ 2004a). Therefore, val- 
ues vary among vegetation types, study sites, and duration times in 
different studies. Furthermore, recent studies have suggested the 


importance of taking both sunlit and shaded foliage into account for 
accurate LUE estimates (He et al., 2013; Middleton et al., 2009) since 
the shaded segment of the canopy is likely to experience less envi- 
ronmental stress (e.g., less solar radiation and lower temperature), 
and hence, has higher LUE. Therefore, an algorithm that integrates 
sunlit and shaded foliage might help to derive more accurate LUE 
values and improve GPP estimates (He et al., 2013). 

5. Summary 

We utilized a simple model with different LUE and fPAR esti- 
mates to simulate GPP for crops and a deciduous forest. Flux tower 
based GPP and fAPAR ch i showed similar seasonal dynamics. The 
MODIS MOD15A2 fPAR showed higher values before green-up 
and after senescence, describing a longer duration of the growing 
season than actually occurred. Estimates using fAPARchi consis- 
tently delivered better agreement with in situ GPP than MODIS 
fPAR. When utilizing s values derived using the MODIS MODI 7- 
like algorithm, underestimations of GPP were observed, especially 
at the agricultural sites. The best performance was found using 
site-specific e max to derive LUE coupled with fAPAR ch i. This study 
demonstrates the importance of accurate land surface parameteri- 
zation to our current carbon monitoring capabilities. In the future, 
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more study sites covering different vegetation types should be 
examined. 
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